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Meyer et al\T\ found that free-standing graphene sheets, just one atom thick, display spontaneous 

ripples. The ripples are of order 2-20 A high and 20-200 A wide. The sheets in which they appear are 

only one atom thick, and extend for around 5000 A through vacuum between metal struts that support 

them. Other groups have since created free-standing graphene as w^ellJ3, Isl, |j], and similar ripples have 

been found for graphene on a glass substrate[5]. Here w^e show^ that these ripples can be explained 

as a consequence of adsorbed molecules sitting on random sites. The adsorbates cause the bonds 

- - - between carbon atoms to lengthen slightly. Static buckles then result from a mechanism like the one 

Q^ that leads to buckling of leaves; buckles caused by roughly 20% coverage of adsorbates are consistent 

3^ w^ith experimental observations. We explain w^hy this mechanism is more likely to explain rippled than 

i^ \ thermal fluctuations or the Mermin- Wagner theorem, which have previously been invoked. 

^^ ■ Fasolino, Los, and Katsnelson[y| and Abdepour et al. [7|| studied the possibility that the ripples are due to thermal 

^ [ fluctuations. It seems unlikely to us that this mechanism explains the observations. Fasolino et al. report that the 

^~5 . bending stiffness k of graphene is on the order of 1 eV. According to Abraham and Nelson[^, the root-mean-square 

QQ ■ fluctuations of a thin sheet are 
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"^ ' where T is the temperature (the experiments are performed at room temperature), and L ~5000 A is the length of the 
sheet. For the purpose of making order of magnitude estimates, we replace KR{q) by the constant k , and get a value 
h «200 A, which is 10 to 100 times larger than the experimental observations. A more vexing problem is created by 
the natural oscillation frequencies of the ripples. The kinetic and bending energy of a thin sheet are approximately 
i 
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Ch ' where p is the mass per area. Thus thermally excited ripples of wavenumber k should be oscillating at a frequency uj 
O , given by 



to' = -k\ (3) 



QQ ■ For waves of scale 200 A , the oscillation frequency comes out to be of order 10^° Hz, while for waves of 50 A the 
ff^ i frequency is higher, 10^^ Hz. These time scales are slow compared with electron dynamics, and thermal oscillations 
On ' of graphene could affect transport properties[Q, [3|. However diffraction through such a rapidly oscillating membrane 
CN ' should be sensitive only to the time average of the density-density correlation function, and that should be completely 
[ — . periodic on time scales much longer than 10"^" s, producing sharp Bragg peaks. However, the analysis of Meyer 
C^ ' et al obtains agreement with the broadened peaks they observed by assuming that the graphene has static ripples. 
00 . Thus the assumption that the observed wrinkling in the graphene layers is purely due to thermal fluctuations seems 
^p^ ' unlikely. 

Another explanation of the ripples[l|, @] invokes the Mermin- Wagner theorem[lO, lll|, according to which a truly 
two-dimensional crystal does not exist except at zero temperature. We do not think that the Mermin- Wagner theorem 
. , is actually signiflcant here. What this theorem says is that when a two-dimensional crystal remains completely planar, 
j^ • it will undergo rotations at large distances that cause a breakdown of long-range order. In the case of graphene at 
room temperature, the scale / on which this happens can be estimated from the expression 

/~aexp[GaVfcBT], (4) 

where a is a lattice spacing and G is the shear modulus of graphene. Since the bulk shear modulus p w G/a of 
graphite is 440 GPa[l2|, and the lattice spacing is greater than a > lA, one has I > lO'^" m. So angular rotation of 
the crystal on the scale of microns should be irrelevant. 
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Figure 1: Energy per area of graphene cylinders plotted versus 1/7?^ where R is the radius of the cylinder. Using (1 0) and 
(1/2 ^/S/4:) as primitive vectors, the honeycomb lattice was wrapped around the y axis. Periodic boundary conditions are 
imposed along the axis direction. 




Figure 2: Ripples in graphene produced by edge effects alone in graphene sheet simulated by MEAM, 100 x 100 A. 



Thus, having decided that neither thermal fluctuations nor the Mermin- Wagner theorem explain rippled graphene 
well, we decided to search for a mechanism that could produce ripples in the ground state. 

Experience with mechanics of thin sheets led us to wonder if graphene might be rippling in a fashion similar to leaves 
and torn plastic|l3l. Il4l. Il5|. In leaves and torn plastic, buckling results from a change in metric. For graphene, this 
mechanism can lead to buckling if for some reason the equihbrium distance between carbon atoms changes at some but 
not all parts of the sheet. To investigate this idea, we used the Modified Embedded Atom Method (MEAM) [la, [l3| 
potential, with parameters of Lee and Lee[18|. To find the bending stiffness, we constructed graphene cylinders of 
various radii R and plotted energy per area versus l/R^, as shown in Figure [TJ The energy per area of a thin elastic 
cyHnder is E = ^n/R^ where k is a bending modulus. We extract a modulus of n =1.77eV from these data, which 
is higher than the value of 1.17eV employed by Fasolino et. uLIq] although of the same order of magnitude. To the 
extent the MEAM potential overstates the bending energy of graphene, it will tend to underestimate the scale of 
ripples. 

A first possibility we checked is that under-coordinated carbon atoms at the edge of the sheet might create ripples 
at the edge that propagate all the way in to the interior of the samples. Figure [2] shows a system 100 A x 100 A in 
size. The edges are free, and the system is initially entirely fiat, except that atoms are randomly displaced from initial 
conditions by around 10^^ A. Then the system is allowed to relax toward a minimum energy state through damped 
dynamics. 

Edge effects do indeed create ripples, and they are of the right size, around 30 A in width and 10 A in height. 
However, the ripples occur only at the edge of the system, and the amplitude decays to zero on the scale of 3.2 A 
(Figure [3]). Thus edge effects alone cannot explain the presence of ripples throughout a sample of size 5000 A. 




Distance into sheet (A) 

Figure 3: Slice through system shown in Figure [2] showing that the ripples decay with a characteristic distance of about 3 A 
away from the edge of the graphene sheet. 




Figure 4: Ripples in graphene produced by 20% coverage of OH, 200x200 A system. 



In correspondence with A. Geim, we learned that OH molecules could be expected to be adsorbed on the surface 
of the graphene sheet. Other molecules such as water and hydrogen may also be present [13, [20, UM- Details of the 
molecule are not Hkely to be critical, and we discuss OH in order to have a specific example. The density of the 
adsorbates is not known, so we have treated it as a free parameter, and considered the effects of randomly placing OH 
molecules on graphene. According to Xu et al[22], attaching an OH molecule to a carbon surface has the consequence 
of increasing the length of two adjoining C-C bonds by around 10% (see their Figure 2, LOD+OH). Rather than 
directly simulating the interactions of graphene with OH, we simply chose carbon atoms randomly from the lattice, 
and increased the equilibrium length to two randomly chosen neighbors by this amount. This corresponds to a small- 
scale change in the metric of the sheet. For the bonds we wished to stretch, the MEAM parameter Re was increased 
to 1.48 A from 1.42 A. This value was chosen because if used for the entire crystal, it produces an equihbrium 
lattice parameter close to what Xu et al find for the C-C bonds stretched by OH. The energy associated with the OH 
attachment is .17 eV. 

As shown in Figure lU with a 20% concentration of OH adsorbates, graphene develops ripples whose wavelength 
and amplitude are comparable to those seen in experiment. Peak-to-peak amplitude of ripples is around six times 
greater than the rms amplitude, so by 40% concentration, the ripples have a peak-to-peak amplitude of around 10 
A. The variation of ripple amplitude and wavelength is displayed in Figure [5j It is not possible to deduce the 
OH concentration from these computations because both experimental and computational uncertainties are too large 
at this point to permit it. Whether this mechanism is in fact responsible for the buckles might be determined by 
experiments on adsorbate-free surfaces conducted in high vacuum. 
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Figure 5: Wavelength and amplitude of buckles in 200 x 200 A sheets as a function of OH concentration. The wavelength changes 
rather little with concentration, while rms amplitude increases. Amplitudes of the rippled peaks are around six times larger 
than the rms amplitude; peak-to-peak amplitude at 40% OH concentration is around 10 A. Wavelength and rms amplituded 
were computed for free-standing sheets after excluding 20 A of material at the edge of the sample. Wavelength was computed 
by decomposing the sheet into a series of line scans, taking the one-dimensional Fourier transform of them in turn, finding 
the average wave vector k for each line weighted by the amplitude of the Fourier transform, computing A = 2-K/k and finally 
averaging A over all the line scans. Error bars represent standard errors after averaging over three independent trials per 
concentration. 
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